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ABSTRACT 

We perform smoothed particle hydrodynamics (SPH) simulations for merging binary carbon-oxygen 
(CO) white dwarfs (WDs) with masses of l.I and I.O Mq, until the merger remnant reaches a dynam¬ 
ically steady state. Using these results, we assess whether the binary could induce a thermonuclear 
explosion, and whether the explosion could be observed as a type la supernova (SN la). We investi¬ 
gate three explosion mechanisms: a helium-ignition following the dynamical merger (‘helium-ignited 
violent merger model’), a carbon-ignition (‘carbon-ignited violent merger model’), and an explosion 
following the formation of the Chandrasekhar mass WD (‘Chandrasekhar mass model’). An explosion 
of the helium-ignited violent merger model is possible, while we predict that the resulting SN ejecta 
are highly asymmetric since its companion star is fully intact at the time of the explosion. The carbon- 
ignited violent merger model can also lead to an explosion. However, the envelope of the exploding 
WD spreads out to ~ O.IRq; it is much larger than that inferred for SN 20IIfe (< O.li?©) while much 
smaller than that for SN 2014J (^ l^^©)- For the particular combination of the WD masses studied 
in this work, the Chandrasekhar mass model is not successful to lead to an SN la explosion. Besides 
these assessments, we investigate the evolution of unbound materials ejected through the merging 
process (‘merger ejecta’), assuming a case where the SN la explosion is not triggered by the helium- 
or carbon-ignition during the merger. The merger ejecta interact with the surrounding interstellar 
medium, and form a shell. The shell has a bolometric luminosity of more than 2 x 10^^ erg s“^ 
lasting for ~ 2 X 10'^ yr. If this is the case. Milky Way should harbor about 10 such shells at any 
given time. The detection of the shell(s) therefore can rule out the helium-ignited and carbon-ignited 
violent merger models as major paths to SN la explosions. 

Subject headings: binaries: close — galaxies: evolution — supernovae: general — white dwarfs — 
hydrodynamics 


1. INTRODUCTION 

The type la supernova (SN la) is one of the brightest 
events in the universe, and plays an important role as 
a cosmological distance indicator. It is widely accepted 
that an SN la is a thermonuclear explosion of a carbon- 
oxygen (CO) white dwarf (WD), and that the explosion 
is triggered by interaction between the CO WD and its 
companion star (Nomoto 1984; Hillebrandt & Niemeyer 
2000). However, it is still controversial whether the com¬ 
panion star is a non-degenerate star (single degenerate 
scenario; SD) (Whelan & Iben 1973; Nomoto 1982), or a 
degenerate star (double degenerate scenario; DD) (Iben 
& Tutukov 1984; Webbink 1984). There are other scenar¬ 
ios, for example, core degenerate (CD) scenario in which 
the companion is an asymptotic giant branch core (Kashi 
& Soker 2011), and collisional DD scenario in which two 
CO WDs collide in a dense stellar cluster, or in a multiple 
stellar system (Aznar-Siguan et al. 2013, 2014). 

The SD scenario has been well tested by recent ob¬ 
servations (see Maoz, Mannucci & Nelemans 2013, for 
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a review). There are multiple observational indications, 
some of which are for the SD scenario and the others 
are against it, for different objects. The observational 
studies against the SD scenario include the following: Li 
et al. (2011) have detected no red giant star in the deep 
pre-explosion images of the site of SN 2011fe. Schae¬ 
fer & Pagnotta (2012) have reported that no main se¬ 
quence or red giant stars are observed at the central 
region of an SN la remnant, SNR 0509-67.5*. On the 
other hand, the observations supporting the SD scenario 
include the following: Dilday et al. (2012) have observed 
SN la PTE llkx, and have found the evidence of a strong 
interaction between the SN ejecta and circumstellar mat¬ 
ter (CSM). Here, the CSM is thought to originate from a 
symbiotic nova, which consists of a WD and a red giant 
star.^ 

* However, we should note that the null detection can be ex¬ 
plained by spin-up/spin-down models, where the companion star 
evolved to become a helium WD during the spin-down phase of 
the CO WD before the delayed carbon ignition occurs in the cen¬ 
ter (Di Stefano, Voss Sz Claeys 2011; Justham 2011; Hachisu, Kato, 
Nomoto 2012). 

^ For CSM in the CD scenario, see Soker et al. (2013). 
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As described above, the SD scenario has been directly 
tested by many studies. On the other hand, most of 
the ‘observational’ support for the DD scenario indeed 
come from ‘non-detection’ - it is supported since the 
SD scenario is ruled out for particular objects. ^ It 
is therefore necessary to assess the DD scenario directly, 
based on theoretical predictions of what should be ob¬ 
served if the DD scenario is the case. For this pur¬ 
pose, the following two questions should be answered: 

(1) Whether a CO WD which accretes materials from a 
companion CO WD results in a thermonuclear explosion. 

(2) Whether such an explosion is observed as an SN la. 

Many previous studies have focused on the first ques¬ 
tion, finding many possible paths in which a primary 
CO WD could explode. These models can be gen¬ 
erally divided into the Chandrasekhar mass and sub- 
Chandrasekhar mass models. In the Chandrasekhar 
mass model, the CO WD reaches the central density 
higher than a critical density to ignite explosive car¬ 
bon burning, whereas it is not the case for the sub- 
Chandrasekhar mass model. 

In the Chandrasekhar mass model, two CO WDs 
merge, and the merger remnant evolves hydrostatically 
after the merger toward an explosion. Whether the 
merger remnant explodes as an SN la depends on the 
structure of the merger remnant. This structure has been 
intensively investigated by means of two different numer¬ 
ical schemes: smoothed particle hydrodynamics (SPH) 
simulations (Benz et al. 1990; Guerrero, Garcia-Berror 
& Isern 2004; Yoon et al. 2007; Loren-Aguilar, Isern & 
Garcia-Berro 2009; Raskin et al. 2012) and mesh-based 
hydrodynamics simulations (D’Souza et al. 2006; Motl et 
al. 2007). Among these works, Yoon et al. (2007) have 
found a path to the SN la explosion. Moreover, Zhu et 
al. (2013) and Dan et al. (2014) have performed large pa¬ 
rameter surveys for various binary CO WD parameters, 
and have searched for the systems which can explode as 
an SN la. 

The sub-Chandrasekhar mass model can be subdi¬ 
vided into several categories. First, Pakmor et al. (2010) 
have suggested a ‘carbon-ignited violent merger model’ 
(see also Pakmor et al. 2011, 2012a,b). In this model, 
hotspots appear in the course of the merger of binary 
CO WDs, and generate carbon detonation leading to an 
explosion. Second, Pakmor et al. (2013) have also pro¬ 
pounded a ‘helium-ignited violent merger model’. In this 
model, a helium layer accreted onto a primary CO WD 
from a companion WD rises in temperature. It is sug¬ 
gested that the helium detonation occurs at some point 
in this layer, and then the shock compression triggers the 
carbon detonation inside the primary CO WD, leading to 
an explosion. This is an analog to the double detonation 
model (e.g. Woosley et al. 1980; Nomoto 1980, 1982), ex¬ 
cept for the nature of the helium donor. Finally, Shen et 
al. (2012), Schwab et al. (2012), and Ji et al. (2013) have 
argued that binary CO WDs could explode shortly after 
the merger due to magnetohydrodynamical effects. 

In this paper, we investigate the helium-ignited vio¬ 
lent merger, carbon-ignited violent merger, and Chan¬ 
drasekhar mass models from two points of view: (1) The 
first point is whether these models lead to a successful 

t See, however, the earlier footnote for the SD scenario and Soker 
et al. (2014) for the CD scenario of SN 2011fe. 


ignition to initiate an explosion. For this purpose, we 
perform SPH simulations of a merger of binary CO WDs. 
However, we can not directly follow initiation of its explo¬ 
sion; in order to follow the initiation, we need SPH sim¬ 
ulation with impossibly high space resolution, say 1 cm. 
Instead, we judge success or failure in the explosion from 
the density and temperature obtained through the SPH 
simulation. Especially, to investigate various possibilities 
in the ignition process, we adopt the following strategy: 
Even when we infer the success of a particular mode of 
the ignition, we do not stop our simulation, and the SPH 
simulation results from the subsequent evolution are used 
to test another mode of the ignition. This is because our 
inference for a particular ignition mode is not decisive, 
and because we want to test several models with different 
ignition conditions. As described above, we use tempera¬ 
ture to infer the success or failure. However, temperature 
is vulnerable to random noises in SPH simulations. In 
order to obtain robust inference, we carefully treat tem¬ 
perature in the following two ways. Eirst, we adopt two 
types of temperatures: raw and smoothed temperatures 
(described in detail later). Second, we do not solve nu¬ 
clear reactions, which are sensitive to random noises, and 
could increase temperature in an unstable manner. 

(2) The second point is whether the expected outcome 
of the explosion is consistent with observations of SNe la. 
For this purpose, we adopt several observational indica¬ 
tions. The first test is the progenitor radii of SN 2011fe 
and SN 2014J. The former and latter radii are inferred 
to be less than O.li?© (Nugent et al. 2011; Bloom et al. 
2012; Zheng et al. 2013; Mazzali et al. 2014) and more 
than li ?0 (Goobar et al. 2015), respectively. As another 
test, we discuss ^®Ni distribution and the ejecta geom¬ 
etry, both of which are thought not to be highly asym¬ 
metric for nearby SNe (Maund et al. 2013; Soker et al. 
2014) and for SN remnants. 

We mainly investigate a merger of a binary consist¬ 
ing of two GO WDs with masses of 1.1 and 1.0 Mq. 
For a benchmark, we also follow evolution of a binary 
with masses of 0.9Mq CO WD and O. 6 M 0 CO WD. 
Other combinations are investigated elsewhere (Sato et 
al. 2015). 

In sum, we find that for the particular binary parame¬ 
ters studied in this paper, the helium- and carbon-ignited 
violent merger models lead to an explosion, while the 
Chandrasekhar mass model does not, according to our 
inference. However, as we mentioned above, our infer¬ 
ence is not decisive due to various uncertainties. We 
therefore suggest a way of evaluating our inference and 
constrain the fates of WD mergers, based on the insights 
obtained through the SPH simulations. If the system in¬ 
deed does not immediately lead to an SN la explosion, 
the binary merger should leave its merger WD remnant 
and the unbound materials ejected from the system dur¬ 
ing the dynamical phase of the merger process (here¬ 
after ‘merger ejecta’). The merger ejecta interact with 
interstellar medium (ISM), and form a shell (hereafter 
’merger shell’), which is analogous to formation of an SN 
remnant. If we find these merger remnants and merger 
shells, we can dismiss our inference and alternatively we 
can use these observational counterparts against the vi¬ 
olent merger scenarios leading to an SN la. We discuss 
the detectability of such events, in particular the merger 
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shell. We thereby suggest that it is possible to detect 
such events. 

This discussion has another benefit. The merger shell 
can be detected not only from pre-explosion images of a 
site of an SN la, but also in the post-explosion observa¬ 
tions, if the SN la happens in the Chandrasekhar mass 
model^. The detection of the merger shell can therefore 
directly support the Chandrasekhar mass model as a re¬ 
sult of the WD merger. 

This paper is structured as follows. In section 2, we 
describe the methods of our simulations. In section 3, 
we show the results from our simulations. In section 4, 
we assess whether our binary CO WDs can explode, and 
whether the explosion can be observed as an SN la. In 
section 5, we discuss the detectability of the merger rem¬ 
nant and merger shell. Finally, we summarize our find¬ 
ings in section 6. 


We use an equation of state (EoS) as functions of pi 
and Ui in order to get not only Pi, but also raw temper¬ 
ature Ti and sound speed Cs,i. For the EoS, we adopt 
the Helmholtz EoS (Timmes & Swesty 2000). This EoS 
includes thermal radiations, an ideal gas of ions, an 
electron-positron gas with an arbitrary degree of relativ¬ 
ity and degeneracy. The EoS requires chemical composi¬ 
tions of fluids. We assume that the chemical composition 
is uniformly fixed to 50 percent of carbon, and 50 percent 
of oxygen in the number fraction. The chemical compo¬ 
sition is fixed during the whole simulations, since we do 
not consider nuclear reactions. 

The gravity is calculated as follows. The gravity Qi is 
the sum of the Newtonian gravity on f-particle exerted by 
all the other particles. We introduce so-called Plummer 
softening to the Newtonian gravity. Thus, it is expressed 
as 


2. METHOD 

In this section, we describe the methods of our simu¬ 
lations. In section 2.1, we briefly present schemes in our 
SPH simulations. In section 2.2, we define several quan¬ 
tities used throughout the paper. In section 2.3, we show 
how to set up an initial condition of binary CO WDs. In 
section 2.4, we summarize a set of physical and numer¬ 
ical parameters used in this study. In section 2.5, we 
introduce our computing environment. 

2 .1. SPH simulation 

We solve Lagrangian hydrodynamics equations with 
self gravity by means of SPH simulations. Each SPH 
particle is evolved by the following equations: 


VP, , 

v^= +g^, 

(1) 

Pt 

zi, = (V • vfi , 

(2) 


Pi 


where Vi, Ui, Pi, and pi are the velocity, specific internal 
energy, pressure, and (mass) density of i-particle, respec¬ 
tively, and Qi is gravity exerting on i-particle. The over¬ 
dots indicate the first time derivative. The nabla symbol 
means an operator of {d/dx,d/dy,d/dz). 

We briefly explain our SPH formulations. In our SPH 
simulations, we solve the ‘vanilla ice’ SPH equations. We 
adopt a cubic spline kernel for the SPH kernel interpo¬ 
lation. The SPH kernel is modified in the same way as 
Thomas & Couchman (1992). Similarly to Rosswog et al. 
(2000), we adopt the treatment of time-dependent artih- 
cial viscosity (Morris & Monaghan 1997), combined with 
a recipe which suppresses the viscosity from shear motion 
(Balsara 1995). This is described in detail in appendix A. 
We set the length of the SPH kernel of a particle, such 
that the arithmetic average number of neighbor parti¬ 
cles over all the particles, (n^), is 150. Hereafter, “()” 
indicates an arithmetic average of quantities over all the 
particles. Neighbor particles of z-particle, Ui, are defined 
as particles whose distances from z-particle is less than 
the kernel length of z-particle. 

§ The Chandrasekhar mass model can be successful when binary 
CO WDs consist of the mass combinations different from the one 
studied in this paper (see Sato et al. 2015). 
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where mj and are respectively the mass and position 
vector of j-particle, G is the gravitational constant, and 
e is the gravitational softening, fixed to e = 3 x 10® cm. 
In practice, we calculate the gravity with an octree al¬ 
gorithm (e.g. Barnes & Hut 1986). In such algorithms, 
gravity exerted on z-particle by distant particles is ap¬ 
proximated as its multipole moment. We consider the 
multipole moment up to dipole moment. In order to 
define whether particles are distant or not, we use Multi¬ 
pole Acceptance Condition (MAC). We choose the same 
MAC as introduced by Salmon & Warren (1994). The 
MAC has one accuracy parameter. A, which is the same 
notation as in Nakasato et al. (2012). The accuracy pa¬ 
rameter A has the dimension of mass divided by square 
of length. As A becomes smaller, gravity is calculated 
with a higher degree of accuracy. At the beginning of our 
simulation, we determine A from Qi at the initial time, 
such as 


A = 


(|g» 

G ’ 


( 4 ) 


where a relative error of gravity on a particle is ensured 
to be less than p. Except for t = 0, we adopt A = 
1.7 X 10^^ gcm~^, which corresponds to zy = 0.01. At 
t = 0, we adopt A ^ 1.0 x 10® gcm“®. The reason why A 
at t = 0 is extremely smaller than A at t > 0 is as follows. 
We do not know the gravity Qi sA t = 0. Therefore, we 
do not know what A corresponds to zy = 0.01 at t = 0. In 
order to avoid obtaining the gravity with a relative error 
larger than 0.01, we use sufficiently small A at t = 0. 
As a result of the calculation of the gravity at t = 0, it 
becomes clear that A ~ 1.0 x 10® gcm“^ corresponds to 
zy = 6 X 10“®. In other words, we calculate the gravity 
with a relative error much less than 0.01 at t = 0. 


2 .2. Definitions 

In this section, we define several quantities used 
throughout this paper. 

We define two types of temperatures: raw and 
smoothed temperatures. The raw temperature of a par¬ 
ticle, Tj-^i, is directly obtained from the Helmholtz EoS 
as a function of pi and Ui, as described in section 2.1. 
Using Ty-^i, the smoothed temperature of a particle, Tg^i, 
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is calculated as 

N 

(5) 

Pj 


where 14^ is the SPH kernel. This smoothed temperature 
is similar to “SPH-smoothed temperature” in Dan et al. 
(2014). 

We define a timescale of a nuclear reaction, tnuc, as 

tnnc = —, (6) 

^nuc 

where Cp is specific heat at constant pressure, Cnuc is an 
energy generation rate per unit mass for the nuclear re¬ 
action, and T is either the raw or smoothed tempera¬ 
ture. We indicate tnuc,r and tnuc,s as the timescale of a 
nuclear reaction, when we calculate the timescale with 
the raw and smoothed temperatures, respectively. We 
replace the subscripts “nuc” in equation ( 6 ) with “3a” 
for the triple-alpha reaction, and “cc” for the -|- 
reaction. We calculate an energy generation rate by the 
triple-alpha reaction in the same way as Kippenhahn & 
Weigert (1990): 

£30 = 9 he/ 3 op"^|T 8-3 X exp (-44.027/r8), (7) 

where qhe = 5.09 x 10^^ [erg g“^ s“^], fza = exp(2.76 x 
10-3pi/2r-3/2) 

is the weak electron screening factor 
(Salpeter 1954; Clayton 1968), X 4 is the mass fraction 
of helium, and Tg = T/ 10 ®, which is the same choice as 
Dan et al. (2014). We also define an energy generation 
rate by -I- as 

Ecc = pgcPc^T9 exp(-Q/Tg/^ -k /cc), (8) 

where qc = 4.48 x 10^® [erg mol“^], fee is a screen¬ 
ing factor (Blinnikov & Khokhlov 1987), Atq = 8.54 x 
10262^5/6^,-3/2 g ^ 84.165, Tg = T/(10®K), 

and Tga = 79/(1 -k O.O 67 T 9 ) (Fowler, Caughlan & Zim¬ 
merman 1975 ). We calculate carbon abundance as Yc = 
nc/ipNa) = 0.033 [g“^], where nc is the number density 
of carbon, and is Avogadro constant. 

We define the local dynamical timescale as 

tdyn = (24^Gp)-^/" , (9) 


which is the same choice as that of Nomoto (1982). 

We introduce a shock detector of a particle in order to 
search for shock-heated regions. The shock detector is 
defined as 


D 


S,l 


-h 


— (V • Vi) 


( 10 ) 


where fi is Balsara switch (see appendix A). This indi¬ 
cates how strongly a fluid element is compressed. The 
element is compressed when Dgy > 0 , and is extended 
when Dg i < 0. The critical value between shocked and 
unshocked regions is about unity, but is not severe. Al¬ 
though it is just an indicator, it should be useful to grasp 
where fluid elements are compressed by shock waves. 

We present the definition of merger ejecta as follows. 
We consider a specific orbital energy of a particle at the 
time t as bi = + 0.5|ui(t)|^, where 0 i(t) and Vi{t) 


are a specific potential energy and velocity of the particle 
at the time t. If a particle has bi{f) > 0, it is unbound at 
the time t. A particle which keeps biit') > 0 during t' >t 
is defined as a merger ejecta at the time t. This means 
that a particle unbound temporarily is not counted as a 
merger ejecta. We also show the definition of a terminal 
velocity of a merger ejecta: 

Uej,i(t) = (26,)^^". (11) 

We can interpret the terminal velocity as a velocity of a 
particle at infinity. 


2.3. Initial condition setup 

We set up an initial condition in our SPH simulations 
in three steps. Our setup method is the same as that 
of Dan et al. (2011), unless otherwise noted. In the first 
step, we generate two single CO WDs individually. In 
the second step, we choose particles from each of the 
single CO WDs, and assume that these particles consist 
of pure helium. In the third step, we combine the two 
CO WDs in the same frame, and form binary CO WDs. 
We explain these steps in detail in the following. 

We can divide the first step into five substeps. In the 
first step, we make 1 -dimensional density profile of a 
fully-degenerate CO WD with uniform temperature of 
10® K. In the second substep, we map particles in such a 
way that their mass densities and specific internal ener¬ 
gies are consistent with the single CO WD. In the third 
substep, we relax these particles as follows. We evolve 
these particles for 20 s in the simulation time by means 
of SPH simulation which is different from that described 
in section 2.1 in two points. One is that a specific inter¬ 
nal energy of each particle is fixed. In other words, we 
do not solve equation (2). The other is that each particle 
receives a damping force against its motion. The damp¬ 
ing force on Aparticle is added to the right-hand side of 
equation ( 1 ), and given by 


{Vi 


)damp 


Vj 

"^damp 


( 12 ) 


where we set T^amp — GdampAt, and Gdamp — 128. We 
adopt Gdamp = 128 experientially. Nevertheless, the 
value of Gdamp does not have the impact on the structure 
of the CO WD. The structure of the CO WD in the case 
of Gdamp = 128 is almost the same as that in the case 
of Gdamp = 64. In the fourth substep, we relax these 
particles again by evolving them during 80 s by means of 
SPH simulation described in section 2.1. The purpose of 
this substep is to avoid mixing of the helium and carbon- 
oxygen particles. This substep has not been done in Dan 
et al. (2011). In the fifth substep, we shift the center of 
mass in the positions of the SPH particles to the origin 
of coordinates, and the center of mass in the velocities 
to zero. 

In the second step, we choose the outermost particles 
of each of the single CO WDs, and assume these particles 
consist of pure helium. We call these particles ‘helium 
particles’. We define a helium fraction, /ne, as the ra¬ 
tio of helium mass to CO mass in these particles. We 
note that the EoS of these helium particles is not helium 
one, but carbon-oxygen one, but it would not introduce 
a large error for our purposes. We compare the radius 
of a pure CO WD with that of a CO WD with helium 
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(/He = 1 X 10-3, 3 X 10-4, and 4 x lO-^) by perform¬ 
ing 1-dimensional hydrostatic calculation, and find that 
the difference among their radii is at most 1 percent. 
Proactively, we show that the helium particles are not 
mixed carbon-oxygen ones without merging. Figure 1 in¬ 
dicates the distribution of the helium and carbon-oxygen 
particles. Particularly, as seen in the companion WD 
which does not accrete materials, the helium particles 
keep staying on the surface of the companion WD. 

We can divide the third step into two substeps. In the 
first substep, we put two single CO WDs in the same 
frame. Hereafter, we call the more massive CO WD 
‘the primary’, and the less massive one ‘the companion’. 
We set a separation between the primary and compan¬ 
ion, such that the companion’s Roche-lobe radius is Ciobe 
times larger than the companion’s radius. We estimate 
the Roche-lobe radius from an approximate formula of 
Eggleton (1983). We can express the separation between 
the primary and companion, ag, as 


ao — Ciobe 


0.49g2/3 


-1 


0 6g2/3 log(l + 


(13) 


where Rc is the companion’s radius, and q is the ratio of 
the companion’s mass (Me) to the primary’s mass (Mp). 
We set Ciobe = 2 in this study. Hereafter, we define the 
separation between the primary and companion as the 
distance between their centers of mass. They are on an 
circular orbit around the origin of coordinates. In the 
second substep, we relax the conhguration of particles 
composing the binary CO WDs, and simultaneously de¬ 
cay the orbit of the binary. For the relaxation, we again 
introduce the damping force given by equation (12). We 
decay their orbit at every Atdecay seconds. The extent 
of the decay is given by 

Audecay — Atdecayj 

'T’decay 


where a is the binary separation before the decay. We 
can write Tdecay as 

Tdecay = y:; , i (15) 

<-/decay a/ I^Pc,0 


where pc,o is the overall mass density of the companion 
before this relaxation and orbital decay, and written as 
Pc.o = Mc/(47ri?3/3/ \Ye adopt Atdecay = 1/64 s and 
Cdecay = 0.05. In this relaxation and orbital decay, we 
take a co-rotating frame of reference. Therefore, we add 
centrifugal and Coriolis terms to equation (I), expressed 
as 


field, where the primary and companion are approxi¬ 
mated to be point mass at their center of mass. The 
potential energies are expressed as 


4>app — 


GMp 

l’’-’’pl 


GMc 

\r - rd 


i (w X r) 


2 


(17) 


where rp and rc are the position of the center of mass 
of the primary and companion, respectively. The red 
curves (<i>nuni) indicate the potential energies of the field, 
where the potential energies among particles are calcu¬ 
lated numerically. The potential energies are expressed 
as 

«>app = $(r) - i (w X r)^ (18) 


where <I>(r) is the potential energies among particles. 
The binary is tidally locked. 


2.4. Simulation run 

We simulate mergers of two types of binary CO WDs: 
pairs of 1.\Mq and I.OMq CO WDs, and O.9M0 and 
O.6M0 CO WDs. The latter is used for a benchmark, 
since such a pair of CO WDs has been widely investigated 
in various studies, such as Yoon et al. (2007). 

We perform the simulations for binary CO WDs with 
masses of I.IM0 and I.OM0 as follows. We set a sep¬ 
aration between the primary and companion, such that 
Ocrit = 1.5 X 109 cm. We choose the helium fraction 
as /ue = 4 X 10-9,3 X 10-^, and 1 x lO-^ for both 
the primary and companion. We adopt various mass 
resolutions, where the numbers of SPH particles used 
to resolve O.IM0 are 64fc, 128A:, 256fc, and 512fc (where 
\k = 2^9 = 1024). The total numbers of particles in each 
model are therefore about 1.4, 2.8, 5.5, and 11 x 10®. We 
name these models ‘model a;M’ {x = 1.4,2.8,5.5, and 
11) after the total number of particles used in each run. 
We follow the evolutions of binary CO WDs for 500 s 
of the simulation time, except for model IIM. In all the 
models, they merge after the binary components orbit 
around each other several times. In all models but for 
model IIM, the simulation is followed until the merger 
remnants reach a dynamically steady state. We stop the 
simulation of IIM just after their merger, since the sim¬ 
ulation time is quite long. 

For binary CO WDs with masses of O.9M0 and O.6M0, 
we separate the primary and companion by Ucrit = 
2.45 X 10® cm. The number of SPH particles to resolve 
O.IM0 is 64A:; the total number of particles is about 
980 X 103. We evolve it for 1000 s. At that time, it 
reaches a dynamically steady state. 


(i9)corotate = X (w X r,) - 2a; X Vi, (16) 

where a; is the angular velocity vector, and |a;| = 
[G(Mp -I- Mc)/a3]^/9^ g^Qp ^j^jg process, Atdecay af¬ 
ter their separation is decreased to less than a critical 
separation Ocrit- This is different from the setup method 
of Dan et al. (2011), who stop this process when any 
particle exceeds the Roche lobe. 

In Figure 2, black dots show the potential energies of 
particles in the frame corotating with the binary. This 
figure corresponds to hgure 3 of Dan et al. (2011). The 
blue curves ('happ) indicate the potential energies of the 


2.5. Computing environment 

We use a code called OTOO (Nakasato et al. 2012), 
which stands for ‘OcTree On Opencl’. The OTOO 
code (hereafter, OTOO) supports a variety of astronom¬ 
ical particle simulations, such as A-body and SPH sim¬ 
ulations. OTOO utilizes an octree algorithm (e.g. Barnes 
& Hut 1986) for fast calculations of particle-particle in¬ 
teractions. It is optimized to multi- and many-core ar¬ 
chitectures on shared-memory environment. It is rela¬ 
tively machine-independent, since it is implemented with 
OpenCL. 
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We perform each SPH simulation on a single node of a 
supercomputer HA-PACS at Center for Computational 
Sciences, University of Tsukuba. A single node of HA- 
PACS consists of two CPUs plus four GPUs. The CPUs 
are Intel Sandy Bridge-EP-8, and the GPUs are NVIDIA 
Tesla M2090. In this configuration, OTOO spends 1.7, 3.3 
6.5, and 13 s for every timestep for models 1.4M, 2.8M, 
5.5M, and IIM, respectively. The total wall-clock time 
for these models are about 160, 360, 670, and 1200 hours. 

3. SIMULATION RESULTS 

In this section, we show results of our simulations. In 
section 3.1, we overview the time evolution of the binary 
CO WDs. In section 3.2, we compare our results with 
previous studies. 

3.1. Overview 

In Figure 3, we present the time evolution of the binary 
CO WDs in models 1.4M, 2.8M. 5.5M, and IIM. In the 
top panels, we show the separations of the binaries. In 
the second top panels, we present mass distribution of the 
binaries. For this purpose, we use ‘x percent Lagrangian 
radii’ defined as follows: Fach of them is a radius of 
a sphere which encloses x percent of the total mass of 
the binary CO WDs, and whose center is the center of 
mass of the primary. In the second top panel, we draw 
10, 50, 60, 65, 70, 75, 80, 90, 99, 99.9, 99.99, and 99.999 
percent Lagrangian radii from bottom to top. In the sec¬ 
ond bottom and bottom panels, we depict the maximum 
Tr.i and Tg^i, respectively, in a range of mass density of 
lO’® < Pi/{g cm“^) < lO’®’*'®'®. These maximum Ty-^i and 
Ts.i are, respectively, indicated as 77, max, a and Ts^ma,x,x- 
We sample all the above quantities at every 1 s. 

We can see in the top panels that the separation is 
steeply decreased at t ~ 270, 180, 120, and 100 s for 
models 1.4M, 2.8M, 5.5M, and IIM, respectively. At 
these times, the binaries merge. We call these times 
‘merger times’. The binaries rotate around each other 
at least 5 times before the dynamical merger, since their 
periods are about 22 s. It is difficult to avoid such differ¬ 
ent merger times for different resolutions, since the mass 
transfer from the companion to the primary is quite un¬ 
stable and chaotic. The separation oscillates with the 
periods of the binaries from the initial time to the merger 
time. This is due to non-zero eccentricities of the orbits 
of these binaries. 

In our simulation, the merger time is relatively smaller 
than previous studies, such as those of Dan et al. (2011). 
This may be because we stop the relaxation process of the 
binary when the binary separation becomes smaller than 
a critical separation. We may make the binary separation 
too small. In section 4, we discuss this effect on whether 
the thermonuclear explosion becomes successful or not. 

We follow the time evolution of the inner mass distri¬ 
bution (< 90 percent Lagrangian radii) of model 5.5M 
as an example. This is instructive, since the mass distri¬ 
bution in other models is similar to that in model 5.5M, 
except their merger times. In model 5.5M, the 10 and 50 
percent Lagrangian radii are not changed throughout the 
simulation. These radii represent the primary’s material. 
This means that the merger has little effect on the mass 
distribution of the primary. On the other hand, the 60 
- 90 percent Lagrangian radii are drastically decreased 


around the merger time (t ~ 120 s). These radii con¬ 
tain the companion’s mass. This means that most of the 
companion’s mass is accreted to the primary at once at 
the merger time. During about 50 s after the merger 
time, the mass distribution of the accreted material is 
still changing. During this time, the 60, 65, and 70 per¬ 
cent Lagrangian radii decrease gradually, whereas the 80, 
and 90 percent Lagrangian radii expand rapidly. After 
this time, the materials below the 70 percent Lagrangian 
radii keep constant in their radii, while the materials 
outer than this in the mass coordinate expand slowly. 
We conclude that the mass distribution, in particular the 
interior of the ‘merger remnant’, is dynamically steady 
at the end of our simulation {t = 500 s). 

Below, the evolution of the outer mass distribution, 
above than 90 percent Lagrangian radii., is explained for 
model 5.5M. The 99.99 and 99.999 percent Lagrangian 
radii rapidly increase well in advance of the merger 
time. Similarly, the 99 and 99.9 percent Lagrangian 
radii rapidly increase around the merger time. The evo¬ 
lution of these Lagrangian radii means that a substan¬ 
tial amount of materials is ejected before/at the merging 
process, which surround the binary system. This matter 
can/should affect observations of outcome of the merger, 
e.g., an explosion, as discussed in detail in section 4. 

The highest raw temperature is achieved for 20 s fol¬ 
lowing the merger time in each model. Hereafter, this 
peak of temperature is called ‘first peak’. The tem¬ 
perature at the first peak is 2.4 x 10® K, 3.6 x 10® K, 
3.8 X 10® K, and 3.8 x 10® K for models 1.4M, 2.8M, 
5.5M, and IIM, respectively. This temperature seems to 
converge to 3.8 x 10® K toward the higher resolution. The 
smoothed temperature also reaches a high value around 
the first peak. In fact, the smoothed temperature at 
the first peak is the highest among the values obtained 
for the whole evolution of the system, in each of models 
2.8M, 5.5M, and IIM. It is 1.4 x 10® K, 1.6 x 10® K, 
and 2.1 x 10® K for models 2.8M, 5.5M, and IIM, re¬ 
spectively. This temperature however does not converge 
even with the highest resolution among our runs. More¬ 
over, this smoothed temperature is much lower than the 
raw temperature at the first peak. If the mass resolu¬ 
tion becomes extremely high, the smoothed temperature 
at the first peak may converge, and may become con¬ 
sistent with the raw temperature at the first peak. We 
do not discuss the convergence of the smoothed temper¬ 
ature anymore in this paper, but we keep in mind that 
we may underestimate temperature when we adopt the 
smoothed temperature. 

Both of the raw and smoothed temperatures in each 
model have another peak about 20 s after the first peak. 
Hereafter, this peak is called the ‘second peak’. This 
peak is achieved for materials in the range of density 10^ 
- 10^ ® g cm“^ in all the models. At the second peak, 
model 1.4M achieves the highest smoothed temperature. 
The raw and smoothed temperature at the second peak 
converge to 2.2 x 10® K and 1.5 x 10® K, respectively. 

Table 1 summarizes the properties of the first and sec¬ 
ond peaks. The first peaks of the raw and smoothed 
temperatures are not coincident. Since the time lag be¬ 
tween the first peaks is at most 5 s, being sufficiently 
smaller than the orbital time (22 s), these peaks appear 
at the same merging phases. This is true for the second 
peaks of the raw and smoothed temperatures, except for 
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TABLE 1 

Properties of the first and second peaks: their type of 

TEMPERATURE, TEMPERATURE, DENSITY, AND TIME. 


Model 

Peak 

Type T [10® K] 

p [10® gem ®[ 

Time [s[ 

1.4M 

1st 

T 

^ r,? 

2.5 

1.7 

284 


1st 

T 

s,? 

1.3 

6.4 

279 


2nd 

T 

r,? 

1.9 

13 

344 


2nd 

T 

s.^ 

1.4 

6.4 

346 

2.8M 

1st 

T 

r,? 

3.4 

2.9 

194 


1st 

T 

s,? 

1.5 

3.6 

195 


2nd 

T 

r,? 

2.1 

13 

230 


2nd 

T 

^ s.^ 

1.3 

5.0 

226 

5.5M 

1st 

T 

r,? 

3.8 

2.9 

134 


1st 

T 

s,? 

1.6 

1.9 

133 


2nd 

T 

r,? 

2.1 

19 

171 


2nd 

T 

^ s.^ 

1.5 

9.2 

173 

IIM 

1st 

T 

^ r,? 

3.8 

3.8 

110 


1st 

T 

s,? 

2.1 

3.7 

111 


2nd 

T 

^ r,? 

2.2 

13 

138 


2nd 

T 

^ s,? 

1.4 

10 

153 


model IIM. However, in model IIM, the smoothed tem¬ 
perature at the time of the second peak of the raw tem¬ 
perature gets as high as ~ 1.4 x 10® K, comparable to the 
smoothed temperature at its second peak. Accidentally, 
the former smoothed temperature is slightly smaller than 
the latter smoothed temperature. 

At t = 500 s, the raw and smoothed temperatures con¬ 
verge to 9 X 10® K and 8 x 10® K, respectively, for materi¬ 
als whose density is in the range of 10® - 10^ g cm“®. As 
for the range of density lO"^ - 10^ ® g cm“®, the smoothed 
temperature converges to 6 x 10® K, while the raw tem¬ 
perature does not. This is because the raw temperature 
in model 5.5M jumps up at t = 350 s. This jump-up is 
due to an effect of an artificial viscosity adopted in the 
SPH simulations. Even a slight amount of the artificial 
viscosity can highly increase temperature of a particle in 
a high density region, since temperature is sensitive to an 
internal energy especially in the high density region. This 
effect is not seen in models 1.4M and 2.8M, and there¬ 
fore the raw temperature in these models rather than in 
model 5.5M should be correct at t = 500 s. Therefore, 
the raw temperature is 8 x 10® K at the density of 10^ 
- 10^ ® g cm“® at t = 500 s. When we compare tem¬ 
peratures of materials at the density < 10^ g cm“® and 
> 10^ g cm“®, the former temperature is higher than 
the latter temperature. This is because a lower density 
region is shock-heated more strongly during the merger 
event. 

Below, we discuss the properties of merger ejecta in 
model 5.5M. Figure 4 shows the time evolution of the 
mass and kinetic energy of the merger ejecta. The mass 
is small at the merger time {t = 120 s), being only 1.3 x 
1O“®M0. The ejecta mass then increases rapidly during 
t = 140 - 150 s. Finally, the ejecta mass reaches 3.9 x 
1O“®M0 at t = 500 s. At t = 500 s, the ejecta stops 
growing in mass. The ejecta mass reaches a constant 
value at about t = 500 s, while the kinetic energy does 
so at an earlier time, about t = 200 s. This is because 
the merger ejecta with higher velocities are formed and 
ejected at an earlier time. Eventually, their total kinetic 
energy at infinity reaches ~ 3.2 x lO"*' erg. 


Table 2 shows the tidal mass, and the masses and ki¬ 
netic energies of the total ejecta in models 1.4M, 2.8M, 
and 5.5M. We can not investigate these properties of the 
total ejecta in model IIM, since we do not follow the 
evolution of this model until a dynamical steady state is 
reached. The properties of the ejecta are independent of 
the mass resolution. 

The merger ejecta can be divided into two groups, 
depending on their formation mechanism. The merger 
ejecta in the first group can be seen in Figure 5. They 
are generated from tidal tails of the primary and com¬ 
panion, and so they are called ‘tidal ejecta’. They be¬ 
come unbound, since they receive orbital angular mo¬ 
menta from rotating bar potential formed by the binary 
system. They are formed only before t = 130 s. The 
total mass of them is 8.6 x 1O“®M0. The tidal ejecta 
become unbound not due to a shock wave, since they are 
far away from shocked regions. The shocked region is 
identical to that creating hotspots at Pakmor’s time (see 
section 4.2). They are not strong shock, since their shock 
detectors are a bit smaller than unity. 

In Figure 6, we show the moment when a part of the 
merger ejecta in the second group are formed. They are 
generated in a shocked region (see the top panels). Here¬ 
after, they are called ‘shocked ejecta’. The shocked re¬ 
gion arises from a collision between the main body of the 
system and a tidal tail. The merger ejecta are prevented 
from traveling toward the directions of the orbital plane, 
which can be seen in the bottom right panel. This is 
because the tidal tail acts as an obstacle. The shocked 
ejecta are formed several times through the above mech¬ 
anism after t = 130 s. They dominate the total mass of 
the merger ejecta. 

The terminal velocities of the tidal ejecta are typically 
3 - 4x 10® cm s“®, similar to the relative velocity between 
the primary and companion just before their merger. 
Those of the shocked ejecta range from 10^ cm s“^ to 
10® cm s“^. These velocities are the largest just after 
the merger, and decrease gradually. 

3.2. Comparison with previous studies 

In this section, we compare our results with those ob¬ 
tained by previous studies. First, we use the results of a 
binary with masses of I.IM0 and I.OMq, and next use 
those of a binary with masses of O.9M0 and O.6M0. 

For a binary with masses of I.IM0 and I.OM0, we fo¬ 
cus on the raw temperature and the mass of the merger 
ejecta at its final state. We also check the smoothed 
temperature, if it is available for comparison in any of 
the previous studies. Not all the previous studies treat 
binary CO WDs exactly with masses of 1.1 and I.OM0, 
and for such cases we adopt, for comparison to work out, 
previous studies with mass combinations similar to the 
one in our simulations. In our results, the maximum raw 
and smoothed temperatures are, respectively, 9 x 10® K 
and 8 X 10® K (see section 3.1). The corresponding tem¬ 
perature in previous studies is as follows: In Dan et al. 
(2014), the maximum raw and smoothed temperatures 
are, respectively, 11 x 10® K and 8.5 x 10® K in the case 
of binary CO WDs with masses of 1.05 and IM0 (see 
their table Al). In Zhu et al. (2013), their maximum 
raw temperature is 9.4 x 10® K in binary CO WDs with 
masses of 1.0 and I.OM0 (see the column of (T’^ax)8 in 




TABLE 2 

The tidal mass, and the mass and kinetic energy of the total ejecta in each model. 


Model 

tidal ejecta mass 

Total mass [Mq] 

Total kinetic energy [erg] 

1.4M 

2.1 X 10~^ 

5.0 X 10-5 

3.0 X 10"‘'' 

2.8M 

2.1 X 10-5 

3.9 X 10-5 

3.1 X lO^'^ 

5.5M 

1.3 X 10-5 

3.9 X 10-5 

3.2 X lO^'^ 

IIM 

- 

- 

- 


their table 2). In Raskin et al. (2012), the maximum 
raw temperature is about 12.5 x 10® K for a binary with 
masses of 1.06 and 0.9QMq (see their figure 9). Our re¬ 
sults are in good agreement with theirs. 

The total mass of the merger ejecta in our simulation is 
3.9 X 10-®Mq. In Dan et al. (2014), it is 7.64 x 
for a binary with masses of I.OMq and I.OMq, 8.50 x 
for 1 .O 5 M 0 and I.OMq, and 1.453 x IQ-^Mq 
for 1.05Mq and 1 .O 5 M 0 . Our mass is consistent with 
those in Dan et al. (2014) to the first order. 

In our results, the total mass of the tidal ejecta is 8.6 x 
10~^Mq. In Raskin & Kasen (2013), the mass of the 
tidal ejecta is 4.7 x 1 O“®M 0 for a binary with masses 
of 1.06Mq and l.OOM©, and 3.3 x 10-®Mq for 1.20Mq 
and I.O 6 M 0 (see their table 1). According to a fitting 
formula of Dan et al. (2014) (see their eq. All), the 
mass of the tidal ejecta is ^ 1 x 1 O“®M 0 in the case of 
a binary with masses of I.IM 0 and I.OM 0 . Our mass is 
smaller than those in these previous studies by an order 
of magnitude. The small merger time of model 5.5M 
should not affect the small mass of the tidal ejecta. In 
model 1.4M, the CO WDs orbit around each other more 
than ten times, but the mass of the tidal ejecta is less 
than 2.1 x 10~^Mq, which is almost the same as in model 
5.5M (see Table 2). Since the tidal ejecta are a minor 
component, the formation of them would be sensitive to 
the detail of the setup of simulations. We do not discuss 
this discrepancy anymore, but we should keep in mind 
that we may underestimate the total mass of the tidal 
ejecta. 

For a binary with masses of 0.9Mq and O. 6 M 0 , we 
show the time evolution of the raw temperature in Fig¬ 
ure 7. The maximum raw temperature after the merger 
event {t = 1000 s) is 6 x 10® K or less depending on the 
position of the materials within the final merger rem¬ 
nant. This is consistent with those found in previous 
studies (Yoon et al. 2007; Zhu et al. 2013; Dan et al. 
2014). In Nomoto, Kamiya & Nakasato (2013), the cor¬ 
responding temperature was reported to be 8.5 x 10® K. 
This is a bit higher than our results although Nomoto, 
Kamiya & Nakasato (2013) adopted the same simulation 
code as ours. We find that this difference comes from 
the difference in the recipes making the initial condi¬ 
tions. We relax a configuration of the binary CO WDs, 
while Nomoto, Kamiya & Nakasato (2013) did not. This 
is consistent with the argument by Dan et al. (2011), in 
which the maximum raw temperature at the final state 
becomes higher without the relaxation. 

For this set of the binary parameters {0.9Mq and 
O. 6 M 0 ), the total mass in the merger ejecta is about 
8.7 X 1 O“^M 0 , which consists mainly of tidal ejecta. The 
ejecta mass of a binary with masses of 0.9Mq and O. 6 M 0 
is much smaller than that of a binary with masses of 


I.IM 0 and I.OM 0 , since the former merger is less vi¬ 
olent than the latter. The corresponding masses are 
2.0 X 1 O“®M 0 in binary CO WDs with masses of 0.96 and 
O.64M0 in Raskin & Kasen (2013), and 1.0 x 1 O“®M 0 
for a combination of 0.9 and 0.65Mq in Dan et al. (2014) 
from their table Al. Our result is almost consistent with 
the previous studies. 

In summary, our raw temperature in the final state is 
consistent with those of the previous studies. The to¬ 
tal mass of the merger ejecta is also in good agreement 
with the previous studies. However, the mass of the tidal 
ejecta is different from those found in previous studies by 
an order of magnitude. Since the tidal ejecta are a minor 
component, the formation of them would be sensitive to 
the detail of the setup of simulations. 

4. ASSESSMENT OE EXPLOSION MODELS 

In this section, we assess explosion models based on 
the results of our SPH simulations. These models are 
the helium-ignited violent merger, carbon-ignited violent 
merger, Chandrasekhar mass (see section 1). We also 
discuss other branches of models which are categorized 
in none of these models. We investigate these explo¬ 
sion models in chronological order: The helium-ignited 
violent merger model (section 4.1), the carbon-ignited 
violent merger model (section 4.2), the other models 
(section 4.3), and the Chandrasekhar mass model (sec¬ 
tion 4.4). In section 4.1, 4.2, and 4.4, we investigate 
whether binary CO WDs can lead to an explosion, and 
whether the explosion can be observed as an SN la. In 
section 4.3, we only discuss what an explosion should 
look like, assuming the explosion is triggered by these 
models. 

4.1. Helium-ignited violent merger model 

In this model, the helium detonation at the surface of 
the primary induces compression of the core materials 
of the primary. This triggers the carbon detonation in 
the core, subsequently leading to an explosion. In order 
to assess a possibility of such a mode in the explosion, 
we provide two check points; whether the helium deto¬ 
nation is initiated, and whether it triggers the carbon 
detonation. 

We define conditions of initiation of the helium detona¬ 
tion. The helium detonation is a supersonic flame pow¬ 
ered by helium burning, i.e. the triple-alpha reaction. In 
order for the triple-alpha reaction to power the flame, the 
triple-alpha reaction needs to proceed faster than cool¬ 
ing due to an adiabatic expansion. In other words, the 
timescale of the triple-alpha reaction should be shorter 
than the local dynamical timescale of the fluid under con¬ 
sideration. Therefore, we adopt t^a < idyn as a criterion 
for the initiation condition of the helium detonation. We 
should keep in mind that this is only a necessary condi- 
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TABLE 3 

The initiation time of the helium detonation, and expected 
®®N l MASS AS A RESULT OF THE DETONATION. 


temperature 

/hb 

time [s] 

5®Ni [Mq] 

raw 

1 X 10"® 

29 

1.1 X 10-5 

raw 

3 X lO-"' 

48 

1.1 X 10-® 

raw 

4 X 10-5 

76 

3.8 X 10-'^ 

smoothed 

1 X 10-5 

64 

1.1 X 10-5 

smoothed 

3 X 10-4 

73 

3.4 X 10-® 

smoothed 

4 X 10-5 

87 

3.8 X 10-'^ 


tion. For a necessary and sufficient condition, a condition 
that the flame propagates supersonically must be satis¬ 
fied. However, we do not consider such a condition in 
this paper. 

Figure 8 shows the time evolution of the total mass of 
helium particles with < ^dyn in model 5.5M. The 
total mass becomes non-zero before the merger time 
{t = 120 s), regardless of a choice of or tza,s, and 
regardless of the value of /hb- Therefore, the system 
generally satisfies the necessary condition to initiate the 
helium detonation before the merger time. 

We define the initiation time of the helium detonation 
as the time after which the total mass of helium particle 
with tza < tdyn raises more than 1Q~^Mq and keeps to 
be so, in order to exclude numerical noises. The initia¬ 
tion time in all the cases is summarized in Table 3. The 
initiation time is earlier when we adopt tza,T: and a larger 
value of /hb- We can also see Figure 9 that a few helium 
particles have tza < ^dyn at the initiation times of the 
helium detonation in all the cases. 

We ignore the contribution of the alpha process, such 
as ^^C(q;, 7 )^® 0 , to the initiation of the helium detona¬ 
tion. This is because the triple-alpha reaction is much 
faster than the alpha process when temperature is less 
than 10® K (e.g. Shen & Moore 2014). The alpha pro¬ 
cess is important for the propagation of the helium deto¬ 
nation, but unimportant for the initiation of the helium 
detonation. 

The small merger time in our simulation may artifi¬ 
cially increase the temperature of the helium particles, 
since the merging process may start suddenly. We thus 
note that the helium detonation might easily be initiated 
artificially. However, we do not consider this effect in the 
following discussion. 

Next, we consider whether the helium detonation can 
trigger the carbon detonation in the core of the primary. 
The helium detonation is thought to trigger the carbon 
detonation in either of two ways. In one way, the helium 
detonation hits the core of the primary, and directly ig¬ 
nites the carbon detonation there. This is called a ‘direct 
drive’ mechanism. In the other way, the shock wave cre¬ 
ated by the helium detonation propagates into the core 
with little or no carbon burning. If the helium detona¬ 
tion region encloses the core, the shock wave then could 
converge somewhere in the core, and could become suf¬ 
ficiently strong to initiate the carbon detonation. This 
is called a ‘converging shock’ mechanism. For the direct 
drive mechanism, a substantial amount of helium is re¬ 
quired (e.g. Moll & Woosley 2013). It is unlikely that this 
mechanism triggers the carbon detonation in our models, 
since /hb is small. We therefore focus on the converging 


shock mechanism, to assess if our system could lead to 
an explosion by this mode, i.e., the helium-ignited violent 
merger. 

For the converging shock mechanism to successfully 
operate, a helium layer needs to enclose the core of the 
primary. We investigate the distribution of helium at 
the initiation time of the helium detonation, as shown 
in Figure 10. In the cases of /hb = 1 x 10“^ and 3 x 
10“^, the helium particles enclose the core of the primary, 
regardless of a choice of the raw or smoothed temperature 
to define the initiation time of the helium detonation. 
However, in the cases of /hb = 4 x 10“®, the helium 
particles are sparse on the orbital plane, again regardless 
of a choice of the raw or smoothed temperature. 

A reason for this sparseness can be explained as follows. 
The mass accretion from the companion to the primary 
proceeds in the following way; the helium particles first 
fall, and subsequently the carbon-oxygen particles do so. 
The carbon-oxygen particles have already started falling 
onto the surface of the primary at the initiation time of 
the helium detonation. These carbon-oxygen particles 
hitting the surface of the primary on the orbital plane 
push the helium particles away. When /hb is smaller, 
the helium particles are distributed more sparsely on the 
orbital plane at the initiation time of the helium detona¬ 
tion for two reasons: First, because of the smaller /hb, 
all the helium particles are more easily pushed away from 
the orbital plane. Second, as /hb becomes smaller, the 
initiation time becomes closer to the merger time (see 
Table 3). As the initiation time is closer to the merger 
time, the carbon-oxygen particles are accreted by the pri¬ 
mary more violently (see Figure 3), pushing the helium 
particles away more easily. 

From the above, we conclude that only in the case of 
/hb ^ 3 X 10“"^, the helium detonation can potentially 
succeed in triggering the core carbon detonation, and 
can lead to an explosion. In Pakmor et al. (2013), the 
helium detonation encloses the core of their primary in 
/hb = 0.01. This result is consistent with our results. 

As described above, our binary can potentially ex¬ 
plode through the helium-ignited violent merger model, 
if /hb ^ 3 X 10“^. Hereafter we consider expected obser¬ 
vational outcome, focusing on brightness just after the 
explosion, hereafter called “early brightness”. In this 
phase, the optical photons are basically powered by the 
thermal energy content produced by the shock heating 
following the SN explosion, in the envelope of the ex¬ 
ploding star. This energy reservoir is sensitive to the 
structure of the binary, especially to the size of the en¬ 
velope of the exploding progenitor. Figure 11 shows the 
matter distribution of our binary system at the initiation 
time of the helium detonation. It has a tidal-tail struc¬ 
ture at t = 29 and 48 s, and a disk structure at t = 64 
and 73 s. The materials spread out to ~ O.li?© away 
from the center of the primary. 

The early brightness is also affected by a radius of 
the companion, through the interaction between the SN 
ejecta and the companion. The radius of the companion 
at the initiation time of the helium detonation is similar 
to that at the initial time. Then, the companion radius 
is about 5 x 10® cm (7 x lO"®!?©). 

Moreover, the early brightness can be partly powered 
by radioactive decay of ®®Ni synthesized by the helium 
detonation. We estimate an amount of ®®Ni as follows. 
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The helium detonation synthesizes ^®Ni when the density 
is higher than a critical density, ~ 10 ^ g cm“^, after the 
helium detonation passes. Taking into account the shock 
compression, we assume that the helium particles with 
Pi > 10 ® g cm“® at the initiation time of the helium 
detonation are converted to ®®Ni. We show the amount 
of ®®Ni as estimated in this way, in all the cases in the 
fourth column of Table 3. 

We qualitatively compare the expected early bright¬ 
ness resulting from the helium-ignited violent merger 
model of our binary with the early brightness of 
SN 2011fe and SN 2014J, in terms of radii of the pri¬ 
mary and companion at the time of the helium ignition, 
and ®®Ni mass after the helium detonation. The expected 
light curve will be much fainter than the light curve of 
SN 2014J. Following the analysis of its early brightness 
(Goobar et al. 2015), SN 2014J is suggested to have ei¬ 
ther a large primary radius (~ li?©), a larger companion 
radius (~ 4/?©), or a large amount of ®®Ni mass near the 
surface (10“®M©). All of these are larger than found in 
the results of our simulation. Another test is provided 
by SN 201 Ife. The expected early brightness may likely 
be brighter than the light curve of SN 201 Ife; the pro¬ 
genitor of SN 2011fe is suggested to have a small radius, 
< O.li? 0 (Nugent et al. 2011; Bloom et al. 2012; Zheng 
et al. 2013; Mazzali et al. 2014). However, we should 
keep in mind that our ‘envelope’ has tidal-tail and disk¬ 
like structures, which has not been taken into account 
in models to connect the size of the progenitor and the 
early brightness. This structure should affect the ex¬ 
pected early brightness; the light curve may be fainter 
than we expect above because of the small opening angle 
of the envelope from the primary, and may be consistent 
with that of SN 2011fe. In order to quantitatively com¬ 
pare the early brightness resulting from our binary with 
those of SNe 2014J and 2011fe, we have to follow the ex¬ 
plosion of our binary by means of numerical simulations. 
This is our future work. 

Another test is provided by properties of SN remnants. 
Papish et al. (2014) have investigated expected proper¬ 
ties of SN ejecta from an exploding CO WD with a he¬ 
lium WD companion. In their model, the separation is 
O.OSi?©, and the helium WD has a radius of 0.02i?©. The 
resulting SN remnant does not have a spherically sym¬ 
metric shape, since the ejecta of the SN prevents from 
moving beyond the helium WD in this direction. We an¬ 
ticipate that an SN remnant resulting from the explosion 
of our binary through the helium-ignited violent merger 
model has a similar shape to that of Papish et al. (2014). 
This is because their and our binaries have similar open¬ 
ing angles of the companion from the primary, i.e., 14 
degree in Papish et al. (2014) and 19 degree in our case. 
The SN remnant with a non-spherically symmetric shape 
is generally not consistent with a large fraction of SN la 
remnants, which tend to be spherically symmetric. 

This section is summarized as follows: When /hb ^ 
3 X 10“"^, our binary can explode through the helium- 
ignited violent merger model. However, the expected 
explosion has different features from SN la in several 
respects. The explosion has the early brightness much 
brighter than that of SN 2011fe, and much fainter than 
that of SN 2014J. Note that the light curve may be 
fainter than we expect, and may be consistent with that 


of SN 2011fe. In addition, the explosion will lead to 
an SN remnant with a non-spherically symmetric shape, 
which is inconsistent with shapes of a large fraction of 
SN la remnants. We note that it is unclear whether the 
other mass combinations can explode through this mode, 
and whether the explosions can be observed as an SN la. 

4.2. Carbon-ignited violent merger model 

If the system survives with no or insufhcient energy 
injection from the helium detonation, there is a chance 
that the carbon detonation initiated at a hotspot leads 
to an explosion, i.e., the carbon-ignited violent merger 
model. A successful explosion in this mode depends on 
whether the hotspot appears in which the nuclear reac¬ 
tion proceeds rapidly to lead to the carbon detonation. A 
necessary condition for this mode to lead to a successful 
explosion is the following: For the carbon detonation to 
take place, the -I- reaction should proceed rapidly 
to lead to thermonuclear runaway. Therefore, a heating 
rate by the -I- reaction should exceed a cool¬ 
ing rate due to an adiabatic expansion. In other words, 
timescale of the -I- reaction should be shorter 
than local dynamical timescale, i.e. tec < Idyn- We call 
materials which satisfy this condition the hotspots. 

We search for hotspots in our simulation at the time 
when any particle has the density of > 2 x 10 ® g cm“® 
and the raw temperature of > 2.5 x 10® K, and at the 
time of the first peak (see section 3.1). The former time 
is the same as the time of the creation of the hotspots de¬ 
fined by Pakmor et al. (2012a), which is based on results 
of Seitenzahl et al. (2009). We call this time “Pakmor’s 
time”. In Figure 12, we show the densities and temper¬ 
atures of particles at Pakmor’s time (left) and the time 
of the first peak (right) in model IIM. Curves in the top 
and bottom panels indicate contours of Ac/idyn- Note 
that the criterion of Pakmor et al. (2012a) are stronger 
than even tec.r < O.ltdyn, since it considers the decay 
of the carbon detonation. Although we use tec.r < Idyn 
for the criterion of the hotspots below, we also consider 
whether the hotspots appear if we choose the constraints 
of Seitenzahl et al. (2009). 

At Pakmor’s time, there are particles satisfying Gc.r < 
^dyn (the top left panel), but none with Gc.s > Idyn 
(the bottom left panel), where tec.r and tcc,s are evalu¬ 
ated using the raw and smoothed temperatures, respec¬ 
tively. On the other hand, at the time of the first peak, 
tec.r < tdyn (the top right panel) and Gc.s < tdyn (the bot¬ 
tom right panel). Therefore, the creation of the hotspots 
at Pakmor’s time depends on the numerical treatment of 
the temperature, while at the time of the first peak the 
hotspots are created robustly irrespective of the treat¬ 
ment of temperature in an SPH simulation. We note 
that, if we use the criterion of Pakmor et al. (2012a) for 
the creation of the hotspots, the hotspots are not cre¬ 
ated even at the first peak from the point of view of the 
smoothed temperature. 

In our simulation, the merger time is relatively shorter 
than previous studies. However, the time to the merger 
does not affect the peak temperature throughout the 
merging process (Pakmor et al. 2012b). Therefore, the 
above discussion whether the hotspots are created is ro¬ 
bust. 

Hereafter we adopt model 5.5M to investigate the for¬ 
mation mechanism of these hotspots. In model 5.5M, 
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Pakmor’s time is t = 118 s as shown in Figure 13. Black 
dots indicate the hotspots. The hotspots are formed as 
follows. Just before t = 118 s, the companion is tidally 
disrupted. Subsequently, a large amount of the disrupted 
debris is rapidly accreted onto the primary. This can be 
seen in the top panel. The center of the primary is lo¬ 
cated on the red region, and the debris extends in the 
direction of the top left from the primary in this figure. 
The accretion of the debris forms a shocked region at the 
surface of the primary. The shocked region can be seen 
from (a;/10®cm, j//10®cm) = (—0.5,—0.5) to (0,0.5) in 
the bottom panel. In this region, the materials are com¬ 
pressed, and the resulting high raw temperatures lead to 
the creation of the hotspots. 

This behavior is qualitatively in agreement with the 
results by Pakmor et al. (2012a) where a system of bi¬ 
nary CO WDs with masses of 1.1 and 0.9Mq, similar to 
our case, is considered. As seen in their figure 1, their 
hotspots are generated at the primary’s surface at the 
time just after the companion is tidally disrupted. 

In order to investigate the effects of the alpha process, 
we draw the positions of the helium particles when /ue = 
4 X 10“®. The reason why we choose /ne = 4 x I0“® 
is that the binary explodes in the case of /ne ^ 3 x 
10“^ before this time through the helium-ignited violent 
merger mode (see section 4.1). As seen in the top panel of 
Figure 13, the helium particles are far away from the hot 
particles indicated by the black dots. The alpha process 
does not affect the carbon-ignited violent merger mode 
at this time. 

Next, we investigate the creation of the hotspots at 
the time of the first peak. Actually, at the time of the 
first peak in model 5.5M, even the particles with the 
highest temperature does not satisfy the necessary con¬ 
dition to form the hotspots if we adopt the smoothed 
temperature. Despite of the absence of the hotspots in 
model 5.5M (for the smoothed temperature), we adopt 
this model for further investigation for the following rea¬ 
sons: We unfortunately do not follow the evolution of 
model IIM until its merger remnant reaches a dynam¬ 
ically steady state, since the simulation of model IIM 
is highly time-consuming. As a result, we can not as¬ 
sess the Chandrasekhar model and investigate merger 
ejecta, using model IIM. In order to assess all the explo¬ 
sion models with the same simulation model, we assess 
the carbon-ignited violent merger model, using model 
5.5M. Indeed, the particles in model 5.5M obtain high 
smoothed temperatures through the same mechanism as 
the hotspots found in model IIM. Regarding these parti¬ 
cles with the high smoothed temperature as the hotspots 
(which should satisfy the condition for the hotspots in the 
corresponding higher-resolution simulation), we investi¬ 
gate the formation mechanism of these particles. 

Figure 14 shows the states of particles around at the 
time of the first peak. The black dots indicate particles 
with the smoothed temperature exceeding 1.5 x 10® K at 
t = 133 s. The number of these particles is 4. We regard 
these particles at t = 133 s as the hotspots. Note that 
these particles are different from those regarded as the 
hotspots at Pakmor’s time. 

Similarly to Pakmor’s time (see Figure 13), we draw 
the helium particles (/ne = 4 x 10“®) on the top right 
panel of Figure 14. The reason why we choose /ne = 


4 X 10“® is the same as the case of Pakmor’s time. Simi¬ 
larly to Pakmor’s time, the helium particles are far away 
from the hot particles indicated by the black dots. There¬ 
fore, the alpha process does not affect the carbon-ignited 
violent merger mode at this time. 

We follow trajectories of these particles. At t = 128 s, 
they are caught between the primary and a tidal tail. 
A collision between the primary and the tidal tail forms 
a shocked region. When the particles pass across the 
shocked region, their kinetic energies are converted to the 
internal energies. After t = 128 s, these particles orbit 
around the primary. At t = 133 s, these particles irrupt 
into a clump created by a debris of the tidally disrupted 
companion. The clump has a high density, 10® ® < 
Pi/(g cm“®) < 10® ®”^®, despite that it is separated from 
the center of the primary by ^ 0.5 x 10® cm. Its surround¬ 
ings have a lower density, 10® < Pi/{g cm“®) < 10® ®. 
When the particles irrupt into the clump, they are com¬ 
pressed (nearly adiabatically), and achieve the highest 
smoothed temperatures. 

Figure 15 shows the creation and evolution of the 
clump quantitatively. The horizontal axes in all but the 
top right panel indicate (p, which is an angle between a 
line segment connecting the coordinate origin and a given 
point, and one connecting the coordinate origin and the 
initial position of the center of mass of the primary. At 
t = 0 s, the densities are almost independent of (p at den¬ 
sity exceeding 10® g cm“®. However, it is not the case at 
t = 130 s, 150 s, and 170 s, during several tens of seconds 
after the merger time. The clump is present at ^ = 5 
(t = 130 s), 4 (t = 150 s), and 0 radian (t = 170 s). At 
t = 250 s, the density becomes independent of (p again, 
announcing that the clump has disappeared. 

Kashyap et al. (2015) have found a hotspot formed 
through a spiral mode instability in the accretion disk 
consisting of the debris of the companion. Such a hotspot 
and spiral possibly appear in our simulation. We can 
see a spiral in the bottom panels of Figure 14. Also, a 
particle with the highest temperature at the second peak 
in our simulation are similar to the hotspot in Kashyap et 
al. (2015); our particle has density of 2 x 10^ g cm“® and 
temperature of 2.1 x 10® K in 5.5M (see Table 1), while 
the hotspot has density of lO"^ g cm“® and temperature of 
3 X 10® K. Our particle has slightly smaller temperature 
than the hotspot, since nuclear reactions are not solved 
in our simulation. We do not discuss the hotspot (or 
the particle with the highest temperature at the second 
peak) anymore. 

Since we consider only the necessary condition to lead 
to the initiation of the carbon detonation and do not deal 
with subsequent evolution following the detonation, it is 
not clear whether this binary system explodes in the end 
or not. Also, it is uncertain whether the explosion is ini¬ 
tiated at Pakmor’s time or at the time of the first peak 
- both are possible but our understanding is currently 
limited by numerical difficulties. Nevertheless, we inves¬ 
tigate what the explosion should look like as we did for 
the helium-ignited violent merger model, especially con¬ 
sidering two situation where the explosion occurs either 
at Pakmor’s time or at the time of the first peak. We 
again focus on the expected early brightness and ®®Ni 
distribution within the hypothesized SN ejecta. 

Figure 16 shows the material distribution of our binary 
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at Pakmor’s time and at the time of the first peak. The 
materials spread out beyond O.li?©, and reach up to ~ 
O. 3 i? 0 . The system has a disk structure, similar to the 
material just before an explosion in the helium-ignited 
violent merger model (see section 4.1). However, the 
disk in this case is more massive and thicker than in the 
case of the helium-ignited violent merger model. 

Following the same argument as presented in sec¬ 
tion 4.1 but applied to the expected pre-SN structure for 
the carbon-ignited violent merger model, we expect that 
this mode results in the early brightness much brighter 
than that of SN 201 Ife, and much fainter than that of 
SN 2014J. This is because the envelope has a radius of 
O. 3 i? 0 , which is much larger than that inferred for 
SN 2011fe (< O.li? 0 ), but much smaller than that for 
SN 2014J (> IRq). The situation is similar to that 
for the helium-ignited violent merger model, and again 
our comparison suffers from the limitation of the spheri¬ 
cally symmetric structure assumed in the estimates of the 
progenitor radii for these SNe. The envelope has disk¬ 
like structure (see Figure 16), and this should be taken 
into account for detailed comparison. Qualitatively, the 
small opening angle of the envelope from the primary 
should make the early brightness fainter than we ex¬ 
pected above, and the light curve may be consistent with 
that of SN 2011fe. This should be quantitatively investi¬ 
gated by means of numerical simulations, similarly to the 
early brightness resulting from the helium-ignited violent 
merger model (see section 4.1). 

Next we discuss ®®Ni distribution synthesized at the 
hypothesized explosion. If the explosion is initiated at 
Pakmor’s time, its distribution is similar to that of Pak- 
mor et al. (2012a). Briefly speaking, ^®Ni is expected to 
be absent in the central region of the SN ejecta in the fol¬ 
lowing reason: Since the companion is burned later than 
the primary, the ashes of the companion are expected to 
be located at the central region in the SN ejecta. At the 
same time, the low density there results in little amount 
of ®®Ni. On the other hand, if the explosion is initiated at 
the time of the first peak, ^®Ni is present at the central 
region of the explosion in the following reason. Since 
the companion has been largely disrupted already, the 
system is more spherically symmetric in this case than 
at Pakmor’s time. The hotspots are created at the sur¬ 
face of the primary, and the explosion is initiated at the 
off-center region. This configuration is similar to that 
in gravitationally confined detonation model (Jordan et 
al. 2008; Meakin et al. 2009) or the off-center delayed- 
detonation model (Kasen et al. 2009; Maeda et al. 2010; 
Seitenzahl et al. 2013), in which a large amount of ®®Ni is 
synthesized near the center of the SN ejecta. The differ¬ 
ence of ®®Ni distribution comes from whether the com¬ 
panion is disrupted or not at the time of the explosion. 

From an observational view point, a model where the 
explosion is initiated at the time of the first peak is more 
favorable, since SN la contains ^®Ni (or other Fe-peak 
elements) near the center of the explosion. However, even 
in this case, after the explosion the ^®Ni distribution is 
expected to evolve to an hourglass-like shape, since ®®Ni 
is prevented from moving toward the direction of the 
orbital plane as blocked by a debris of the companion 
(Raskin et al. 2014). This is not consistent with the 
®®Ni distribution generally inferred for SN la which is a 


spherically symmetric shape (Maund et al. 2013; Soker 
et al. 2014). 

This section is summarized as follows: In the carbon- 
ignited violent merger model, we confirmed that the 
hotspots appear. Therefore, an explosion can poten¬ 
tially occur through this mode. However, it is expected 
that the resulting early brightness is much brighter than 
SN 2011fe and much fainter than SN 2014J. Note that 
the small opening angle of the envelope from the pri¬ 
mary should make the light curve fainter than we expect, 
and that the light curve may be consistent with that of 
SN 201 Ife. Moreover, it is expected that the ®®Ni dis¬ 
tribution does not have a spherically symmetric shape, 
also being inconsistent with ^®Ni distribution generally 
inferred for SN la. The expected observational outcome 
will apply to any explosions through this mode, since the 
explosions can potentially occur only from the WD mass 
combinations similar to one studied in this paper (Sato 
et al. 2015). 

4.3. Other models 

An explosion in the carbon-ignited violent merger 
model happens around the merger time, while an explo¬ 
sion in the Chandrasekhar mass model happens ~ 10"^ 
years after the merger time (e.g. Yoon et al. 2007). Be¬ 
tween these two epochs, an explosion through other mod¬ 
els possibly happens, given that the system does not ex¬ 
perience the explosion through the helium- and carbon- 
ignited violent merger models. For example, Schwab et 
al. (2012) and Ji et al. (2013) have suggested an explo¬ 
sion triggered by magnetohydrodynamical effects in this 
phase. In this section, we do not assess whether the ex¬ 
plosion happens in such a model as it is beyond what we 
can discuss based on our pure hydrodynamic simulation, 
but discuss what the explosion should look like, assuming 
that the explosion does happen. 

The appearance of the explosion will be affected by 
the nature of merger ejecta. Figure 17 shows the spa¬ 
tial and velocity distributions of the merger ejecta at 
500 s. The merger ejecta spread almost isotropically, ex¬ 
cept that they are relatively deficient on the orbital plane 
(see also Figure 6 ). Nevertheless, the covering factor of 
the merger ejecta around the merger remnant is almost 
unity. 

Raskin & Kasen (2013) have also studied effects of 
merger ejecta on observations of a putative SN taking 
place in this phase/mode. However, they have focused 
only on the tidal ejecta for the merger of binary CO WDs 
with masses of O.96M0 and O.64M0. On the other hand, 
our binary model consists of I.IM 0 and I.OM 0 CO WDs. 
Since our binary has more massive components and a 
mass ratio closer to unity, our binary merges more vi¬ 
olently than theirs (Marsh, Nelemans, Steeghs 2004). 
Consequently, our merger ejecta are dominated by the 
shocked ejecta (see section 3.1). Therefore, the expected 
effects of the merger ejecta on the observational features 
are partly different from those described by Raskin & 
Kasen (2013). 

Raskin & Kasen (2013) have shown that NalD absorp¬ 
tion features are potentially observed with a probability 
of 10 ~ 50 percent if the explosion occurs 10® s - 10^ yr 
or 10® yr - 10® yr after the merger time. The probability 
corresponds to the covering factor of their merger ejecta. 
However, in the case of our binary, the NalD absorp- 
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lion features are potentially observed in almost all the 
cases.^ This result is complementary with the discussion 
in section 3.4 and 3.6 of Raskin & Kasen (2013). 

Finally, we point out that, if an explosion does happen 
at t = 500 s, the explosion might look like SN 2014J 
in the early brightness. At t = 500 s, the envelope of 
our binary model spreads out beyond l.Oi?©, as seen in 
Figure 18. The horizontal axis indicates the spherical 
radius. 

4.4. Chandrasekhar mass model 

In this section, we explore a possibility of an explo¬ 
sion in the Chandrasekhar mass model, which might 
take place if the system does not undergo an explosion 
through the mechanisms investigated by the previous sec¬ 
tions. After the remnant reaches a dynamically steady 
state, it gradually loses the thermal energy by neutrino 
cooling, and increases its central density and tempera¬ 
ture. If its central density and temperature exceed criti¬ 
cal values as set by the balance between the nuclear re¬ 
action timescale and dynamical or convection time scale, 
the remnant is likely to explode as an SN la (‘the Chan¬ 
drasekhar mass model’) assuming the central region con¬ 
sists of carbon. On the other hand, it likely collapses 
to a neutron star as triggered by the electron capture, 
if the remnant has become an oxygen-neon-magnesium 
WD before reaching to this phase. 

The presence of carbon in the remnant depends on 
whether carbon is quiescently burned. All of carbon in 
the remnant are burned, if the ^^C -I- ^^C reactions pro¬ 
ceed faster than the neutrino cooling at the time when 
the remnant reaches a dynamically steady state (Saio & 
Nomoto 1985, 1998, 2004). 

Figure 19 shows mass density and temperature of par¬ 
ticles at t = 500 s in model 5.5M. At that time, the 
merger remnant reaches a dynamically steady state. In 
many particles, the timescale of the ^^C -I- ^^C reaction is 
shorter than the timescale of the neutrino cooling, where 
we calculate the neutrino cooling rate using a public code 
available at F. X. Timmes website** which is based on 
Itoh et al. (1996). However, the timescale of -I- 
reaction is longer than the dynamical timescale. There¬ 
fore, the -I- reaction does not trigger the car¬ 
bon detonation, but will convert the merger remnant to 
an oxygen-neon-magnesium WD on thermal timescale. 
We conclude that, after the remnant evolves on thermal 
timescale, the remnant would not explode in the Chan¬ 
drasekhar mass model, rather collapses to a neutron star. 

We should keep in mind that this result may be affected 
by the short merger time in our simulation. According 
to Dan et al. (2011), the maximum temperature in the 
merger remnant becomes high when the binary suddenly 
merges. Therefore, if the time to the merger is longer, the 
maximum temperature in the merger remnant may be 
lower, and the -I- reaction rate is lower than the 
neutrino cooling rate. Then, the remnant could explode 
in the Chandrasekhar mass model. 

5. DETECTABILITY OF MERGER REMNANT AND 
MERGER SHELL 

^ Adding to this, one has to consider the thermal condition of 
the ejecta, which is beyond the scope of this paper. 

** http://cococubed.asu.edu/code_pages/nuloss.shtml 


If binary CO WDs (or its merger remnant) fail to 
explode in all but the Chandrasekhar mass model, the 
merger ejecta have at least 10^ yr of time during which 
the ejecta expand into the ISM. An expanding shell is 
formed as the merger ejecta sweep up its surrounding 
ISM, analogous to an SN remnant. Also, the merger rem¬ 
nant still exists during this phase, analogous to a neutron 
star in a core-collapse SN remnant. For this situation, 
we estimate a detectability of the merger remnant and 
merger shell. In particular, we focus on the merger shell, 
estimating its luminosity. 

The evolution of a merger shell is divided into three 
phases. The first phase is a ‘free expansion phase’, where 
the amount of ISM swept up by the merger ejecta is neg¬ 
ligible. Once the shell sweeps up ISM mass comparable 
to the mass of the merger ejecta, the shell is substantially 
decelerated entering into the ‘Sedov phase’. Materials are 
thermalized behind the shock wave. Because of the high 
temperature, the shell loses only a negligible fraction of 
energy through radiation approximately conserving the 
total energy content. Once the shell is cooled down, a 
signihcant amount of energy is lost by radiation, where 
the momentum is approximately conserved. This phase 
is called a ‘snowplow phase’. 

A luminosity of the shell is notated by Tsheib and given 
by 


Tshell — Akshell, (1^) 


where A is a cooling function, and I4heii is the volume of 
the shell. The volume I4heii can be expressed as Dsheii = 
47ri?g Ashelb where Rs is a radius of the shock wave in 
front of the shell, and Asheii is the thickness of the shell. 
We give the cooling function A as follow: 


A — /^^shell n'e,shelin'H,shell: (^b) 


where rie,shell and nn,shell are the number densities 
of electrons and hydrogen atoms in the shell, Tgheii 
is the temperature in the shell, and (3 = 1.7 x 

10“^® ergcm“^s“^K°-^ (Draine 2011). Note that this 
cooling function takes into account metal emission lines, 
and is applicable in the range of 10® < Tsheii/K < 
10^ ®. Since hydrogen atoms are perfectly ionized after 
they pass through the shock wave, the number densities 
Ue,shell and riH.sheii can be expressed as nH,ism(7 + l)/(7“ 
!)■ 

We estimate the shell luminosity at the free expansion. 
At the free expansion phase, the radius of the shock wave 
Rs oc t, and the thickness Rs is constant. Then, I4heii oc 
t®, and the temperature in the shell, Tgheib is constant. 
Eventually, Tsheii c>c . 

We focus on the Sedov phase, since the luminosity 
Dsheii reaches to the peak around the end of the Se¬ 
dov phase (see below). Although, at the free expansion 
phase, we show only the proportional relation between 
Tsheii and t, we give the equation expressing the relation¬ 
ship between Lgheii and t at the Sedov phase. Assuming 
a strong shock wave, the radius {Rs) and the thickness 
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(Asheii) can be written as 


Solving equation (29), we obtain tcooi as follows: 


- 7 ^ (7-l)(7+l)^ gk.ej 1 ^0.4 

IGtt (37 - 1) mnnisinj 

( 21 ) 


Ashell — 


( 7 - 1 ) 

3(7 + 1) 


Rs 


( 22 ) 


where i?k,ej is the total kinetic energy of the merger 
ejecta, riism is the number density of the ISM in front 
of the shock wave, toh is the mass of a hydrogen atom, 
and 7 is the adiabatic index (Cavaliere & Messina 1976). 
Note that we assume that the ISM consists only of hy¬ 
drogen atom here to provide a first order estimate. 

The shell temperature Tsheii can be obtained as fol¬ 
lows. From Rankine-Hugoniot conditions, the shell pres¬ 
sure psheii is given by 


Pshell — /(7 + 1), (^3) 

where Rg is the speed of the shock wave. On the other 
hand, the pressure Psheii can be written from an equation 
of state of an ideal gas as 


Pshell — ^H,shell^B7shell, 

where fee is the Boltzmann constant. Using equation (23) 
and (24), we obtain the shell temperature: 

T.heii = (25) 

(7 + 1) 

At the Sedov phase, the speed of the shock wave Rs can 
be written as 


R. = 


3 (7-l)(7 + l)^ ^k,e 
47r (3j - 1) TOHUism 


n 0.5 


K 


-1.5 


(26) 


(Cavaliere & Messina 1976). 

Then, we can rewrite equations (19) and (21) as 


L^keii = 2.6 X 10^6 

A-k.ej 


1.68 


3.2 X lO^^erg 


0.32 / ^ \ 2.04 


104yi 




- 0.2 


Icm 

A’kjej 


3.2 X lO'^’^erg J \ lO'^yr 


0.2 / ^ X 0.4 


[ergs ^], 
(27) 


[pc]. (28) 


Cool = 1.1 X 10^ 


0.553 




3.2 X lO'^’^erg) 


0.224 


[yr 


(30) 


The cooling function A is appropriate only when 10® < 
Tsheii/K < 10^®. From equation (25), we obtain the 
following expression: 

T,hen = 2.3xl0®(^^)”°'' 

K 3.2x1ok..g )“(T0^)'”W' m 

From equation (31), the cooling function can be applied 
from a few 10 ^ yrs to a few 10 "^ yrs. 

At the snowplow phase, ISM which passes through the 
shock wave emits energy almost instantly. The shell lu¬ 
minosity is written as Lsheii oc msCs, where rhs is the 
rate of the ISM mass passing the shock wave, and Cg is 
the specific energy which the ISM gains from the shock 
wave. Since Rg oc and Rg oc , rhg oc 
From Rankine-Hugoniot conditions, Cg oc oc 
Then, Lgheii oc at the snowplow phase. 

In Figure 20, we illustrate the time evolution of the 
shell luminosity for a set of typical parameters corre¬ 
sponding to our system; Ak,ej = 3.2 x 10^^ erg, assuming 
iiism = 1 cm“®. We apply equation (27) for the Sedov 
phase. We dehne the time when the Sedov phase be¬ 
gins as the time when the shell has swept up the ISM 
mass comparable to the ejecta mass. The speed of the 
shock wave is set to the average velocity of the merger 
ejecta, 3 x 10® cm s“^ at the free expansion phase. For 
the time when the Sedov phase ends, we adopt tcooi in 
equation (30). 

As seen in Figure 20, the shell luminosity reaches to 
the peak (Tgheii.peak) at the time when the Sedov phase 
ends, tcooi- Substituting tcooi in equation (30) into t in 
equation (27), we obtain the dependence of the peak lu¬ 
minosity of the shell on the kinetic energy of the merger 
ejecta as 

Ashen.peak = 3.1xl0®®(^^)''' 

/ P \ 0 

x(--) [ergs-M. (32) 

V3.2 X 1047erg7 ^ ^ ^ ^ ^ 

We define a lifetime of the shell, Tnfe, during which the 
shell has more than half of the peak luminosity. Then, 
the lifetime is given by 


where we set 7 = 5/3. As seen in the power of t in equa¬ 
tion (27), the shell luminosity keeps increasing as time 
goes by. We estimate roughly the time when the Se¬ 
dov phase is terminated, tcooi, using the shell luminosity 
Lgheii and the kinetic energy of the merger ejecta Ek,ej 
as: 

P^cool 

Ek,ej = / Rshelldt. (29) 

Jo 


Tjife = 4tcool. (33) 

The above estimate on a peak luminosity and lifetime 
of a merger shell is based on results of a merger of WDs 
with masses of \.\Mq and I.OMq. However, such a mas¬ 
sive binary system is rare. To connect the predictions 
to observations, we therefore scale the above result to 
a more common situation of a binary whose total mass 
just exceeds the Chandrasekhar mass. Specifically, we 
consider a merger of WDs with masses of 0.9Mq and 



15 


0.6Mq. This system is representative of the one in which 
the merging process proceeds in the least violent man¬ 
ner amongst the systems potentially leading later to an 
explosion through the Chandrasekhar mass model in the 
ignition mode. According to our simulation for merging 
binary CO WDs with masses of 0.9 and O. 6 M 0 (see sec¬ 
tion 3.2), its merger ejecta have the total kinetic energy of 
7 X lO"*® erg. Consequently, its merger shell will have the 
peak luminosity of Tsheii.peak = 2 x 10 ^® erg s“^, and the 
lifetime of Tufe = 2 x 10'^ yr. We regard these values as the 
lower limits of the peak luminosity and lifetime of merger 
shells, which would be created by any combinations of 
binary WD masses if the total mass exceeds the Chan¬ 
drasekhar mass (hereafter ‘super-Chandrasekhar bina¬ 
ries’). 

The shell could be more luminous than estimated 
above. This is because the shell could be illuminated 
by the merger remnant. Since the merger remnant has a 
high temperature for a while after the merger, it would 
also have a substantial luminosity. This is analogous to 
a planetary nebula. However, we do not consider this 
effect, but simply note that our estimate on the merger 
shell detectability should be regarded as a lower limit. . 

The merger shells will emit ultraviolet and soft X-ray 
photons. We estimate the number of the merger shells in 
Milky Way, iVsheii, assuming that super-Chandrasekhar 
binaries explode only in the Chandrasekhar mass model, 
or fail to explode in all the explosion models. Then, 
When is given by 


When "^10 




M, 


MW 


6 X lOiOM©/ V2 X 104 yr 


Tlife 


(34) 


where Fmerge is a merger rate of super-Chandrasekhar 
binaries per unit mass in Milky Way, Mmw is the mass of 
Milky Way. We adopt Tmerge in Badenes & Maoz (2012), 
and Mmw in McMillan (2C)11) and Licquia & Newman 
(2014). 

We note difficulties in distinguishing merger shells from 
SN remnants and nova shells. Since their explosion en¬ 
ergies are different by several orders of magnitudes, their 
sizes and luminosities can be distinguished, if the number 
density of ISM is known. However, the number density is 
usually unknown. One way to overcome this difficulty to 
identify the merger shells is to use the information about 
the central compact object. 

We search for any hint of the merger shell in the liter¬ 
ature dealing with pre-explosion images of SNe la in or¬ 
der to check whether these SNe la involve merger shells. 
Nielsen et al. (2012) (see also Liu et al. 2012; Nielsen et 
al. 2014) have constrained the upper limit of bolometric 
luminosities of nearby SNe la with pre-explosion images. 
The most stringent limit is on SN 2011fe, ^ 10^® erg s”^. 
This is still lager than the estimated luminosities of 
merger shells. Unfortunately, these images are not so 
deep to constraint the presence of the merger shells. 

As seen in equation (28), a merger shell has a parsec- 
scale size. When an SN occurs within the merger shell, 
the merger shell is not disturbed by the SN ejecta dur¬ 
ing the first 10^ yr after the SN explosion. Therefore, 
existence of the merger shell could be also tested by ob¬ 


servations of an SN la after the explosion. We postpone 
such a study to future. 


6. SUMMARY 

We have performed SPH simulations for merging bi¬ 
nary CO WDs with masses of 1.1 and 1.0 M©, un¬ 
til the merger remnant reaches a dynamically steady 
state. Using these results, we assess whether the binary 
could induce a thermonuclear explosion, and whether 
the explosion could be observed as an SN la. We in¬ 
vestigate three explosion mechanisms: a helium-ignition 
following the dynamical merger (‘helium-ignited violent 
merger model’), a carbon-ignition (‘carbon-ignited vio¬ 
lent merger model’), and an explosion following the for¬ 
mation of the Chandrasekhar mass WD (‘Chandrasekhar 
mass model’). In addition to the evaluation if the result¬ 
ing system satisfies requirements set in each mode. We 
have discussed whether the resulting explosions, through 
different ignition modes, would look like SNe la. 

Our results are summarized as follows: 

• In the helium-ignited violent merger model, our bi¬ 
nary can explode, if the mass fraction of helium ex¬ 
ceeds a critical value, i.e., /ne ^ 3 x 10“4. However, 
the expected early brightness is likely different from 
those of SN 2011fe and SN 2014J, since materials of 
our binary spread out to ^ O.li?©, which does not 
fit to what were inferred for these SNe. Moreover, 
the explosion likely results in an SN remnant with 
an extremely asymmetric symmetric shape, which 
is unusual for SN la. 

• In the carbon-ignited violent merger model, our bi¬ 
nary can explode. However, the explosion likely re¬ 
sults in the early brightness dissimilar to those of 
SN 20IIfe and SN 20I4J for the same reason as for 
the helium-ignited violent merger model. More¬ 
over, we predict that the explosion will synthe¬ 
size ®®Ni whose distribution is extremely aspher- 
ical. This is inconsistent with ®®Ni distribution 
generally inferred for SN la. Note that ^®Ni dis¬ 
tribution will depend on our choice of the raw and 
smoothed temperatures for the ^^C -|- ^^C reaction. 
In the case of the raw temperature, ^®Ni is absent 
in the center of the explosion. On the other hand, 
in the case of the smoothed temperature, ®®Ni dis¬ 
tribution is similar to a hourglass. 

• If our binary explodes a few hundred seconds af¬ 
ter its merger by some mechanism (i.e., the ‘other 
model’), the explosion may have the early bright¬ 
ness consistent with that of SN 20I4J. At that time, 
materials of our binary spread out beyond li?©. 

• For a particular set of the binary parameters ex¬ 
amined in this paper, the binary would not lead 
to an SN la explosion through the Chandrasekhar 
mass model. Rather, the merger remnant should 
be converted to an oxygen-neon-magnesium WD, 
and then will experience an accretion-induced col¬ 
lapse to become a neutron star. 

Binary CO WDs generate the merger ejecta before 
and after its merger time. The merger ejecta will inter¬ 
act with its surrounding ISM, and form a merger shell. 



16 


We estimate a bolometric luminosity of the merger shell; 
the luminosity is more than ~ 2 x 10^® erg s“^ at its 
peak, if the total mass of the binary CO WDs exceeds 
the Chandrasekhar mass. Suppose that all the super- 
Chandrasekhar binaries explode in the Chandrasekhar 
mass model or fail to explode at all, the number of the 
merger shells in Milky Way is estimate to be ~ 10. De¬ 
tection of such merger shells can rule out the helium- 
ignited and carbon-ignited violent merger models. If 
an explosion is initiated in the Chandrasekhar mass 
model, a merger shell can be detected not only from pre¬ 
explosion images of a site of an SN la, but also in the 
post-explosion observations. Unfortunately, we have not 
found merger shells from pre-explosion images of previ¬ 
ous SNe la, since the lower limit in these observations 
is at best ~ 10^® erg s In future, the merger shells 
would be found from post-explosion images of sites of 


nearby SNe la. 
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APPENDIX 

ARTIFICIAL VISCOSITY 

We describe our chosen artificial viscosity, and its numerical parameters. In this section, variables are defined again 
with the same notations as those in the main text. The viscosity exerted on i-particle by j-particle is indicated as H^. 
As described in the main text, we adopt a sheer-free viscosity term (Balsara 1995) combined with time dependent 
viscosity parameters (Morris & Monaghan 1997). Then, the term is given by 


where fij is so-called Balsara switch, and is a sort of bulk and von-Neumann-Richtmyer viscosities. Hereafter, 

we define a combination of the overline and subscript “ij” as 


a:, 




(A2) 


The viscosity term Hy^max is expressed as 





Pij 


{nj ■ V,j < 0) 

(r„ • V,j > 0) 


(A3) 

(A4) 


where pi and Cgy are the mass density and sound speed of i-particle, respectively. The vectors and Vi are, respectively, 
the position and velocity of f-particle, and = rj — ri and Vij = Vj — Vi. The variable pij is given by 


pij — 


hij T ij ' Vij 
|ry |2 -p Sihij ’ 


(A5) 


where hi is the kernel length of f-particle, and = 0.01 is chosen. The viscosity parameter ai is time-dependent. It 
is evolved as 


a 


i — 


^min 

hi/ 


-Pmax[-(V • Ui)(amax 


Vli ), 0] , 


(A6) 


where amax = 1-5, amin = 0.05, and ^ = 0.25. Another viscosity parameter Pi is proportional to such that Pi = 2ai. 
The Balsara switch can be written by fi and fj with equation (A2), and fi is given by 




__ 

|V • Vi\ -P |V X Vi\ -P bic^^i/hi ' 


(A7) 


where we adopt = 0.0001. 
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Fig. 1.— Distribution of helium and carbon-oxygen particles at t = 0, 10, 30, and 50 s in model 5.5M. The red and blue points indicate 
the helium particles of the primary and companion WDs, respectively. The yellow and light blue points show the carbon-oxygen particles 
of the primary and companion WDs, respectively. 

































































































































































































R}tential energy [10 


18 



Fig. 2.— Potential energies of particles as a function of their a:-coordinate. The definition of ^app $num ^-re in the main text. 
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Fig. 3.— Time evolution of binary CO WDs in models 1.4M, 2.8M, 5.5M, and IIM from left to right. Each top panel shows the 
separation between the primary and companion. Each second top panel draws 10, 50, 60, 65, 70, 75, 80, 90, 99, 99.9, 99.99, and 99.999 
percent Lagrangian radii (defined in the main text) from bottom to top. The second bottom and bottom panels indicate the maximum of 
raw and smoothed temperatures in a range of mass density shown in the left panels. 
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Fig. 5.— Shock detector distribution at t = 120 s. Particles with \z\ < 10^ cm are drawn. Black dots indicate merger ejecta at t = 120 s. 
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Fig. 6.— Shock detector distributions at t = 128, 130, and 132 s. Particles with \z\ < 10^ cm and with \y\ < 10^ cm are drawn in the 
top and bottom panels, respectively. Black dots indicate merger ejecta at t = 130 s. 
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Fig. 7.— Time evolution of maximum raw temperature at ranges of mass density in the case of binary CO WDs with 0.9 and O.OMq. 
The ranges of mass density are indicated in the panel. 
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Fig. 8.— Time evolution of the total mass of helium particles with tsa,T < ^dyn (fop) i 3 a,s < ^dyn (bottom). 
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Fig. 9.— Mass density and temperature of helium particles at the time when the helium detonation is initiated. In the left and right 
panels, vertical axes indicate the raw and smoothed temperatures, respectively. From top to bottom, /He = 1 x 10“^, 3 x 10“^, and 
4 X 10“^. A curve in each panel shows mass density and temperature at which timescale of the triple-alpha reaction is equal to local 
dynamical timescale. 
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Fig. 10.— Distribution of helium particles with |y| < 2 x 10® cm at the initiation times of the helium detonation. Red and blue points 
indicate helium particles originating from the primary and companion, respectively. From top to bottom, /He = 1 x 10~®,3 x 10“^^ and 
4 X 10“^. In the left and right panels, the raw and smoothed temperatures are chosen for the initiation conditions of the helium detonation, 
respectively. 
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Fig. 11. — Material distribution of a binary with masses of I.IMq and I.OMq in model 5.5M. We draw only particles separated from 
the center of the primary by 2 x 10^ cm. The time is indicated at the top right in each panel. The time is the initiation time of the 
helium detonation in the cases of /ne = 1 x 10“^ and the raw temperature (two top-left panels), /jje = 1 x 10“^ and the smoothed 
temperature (two top-right panels), /He = 3 x 10”“^ and the raw temperature (two bottom-left panels), and /hc = 3 x 10~^ and the 
smoothed temperature (two bottom-right panels). Dashed curves indicate 0.1i?Q from the center of the primary. 
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Fig. 13.— Mass density (top) and shock detector (bottom) of particles at 118 s in model 5.5M. The coordinate origin is at the center of 
mass of the binary system. The particles are separated from the orbital plane {x-y plane) by < 0.1 x 10^ cm. Black dots indicate particles 
with pi > 2 X 10® g cm“^ and > 2.5 x 10^ K. In the top panel, black crosses indicate the helium particles when /jje = 4 x 10“®. 
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Fig. 14.— Mass densities (top) and shock detectors (bottom) of particles at t = 128 s, 130 s, and 133 s from left to right in model 5.5M. 
The particles are separated from the orbital plane by < 0.1 x 10^ cm. They are colored in the same way as Figure 13. Black dots indicate 
particles with smoothed temperatures more than 1.5 x 10^ K at t = 133 s. The numbers of these particles are 4. In each panel, black 
dots look like only one dot. This is because these particles are located on almost the same positions. In the top right panel, black crosses 
indicate the helium particles when /hc = 4 x 10“ . 
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Fig. 16. — Material distribution of a binary with I.IMq and I.OMq in model 5.5M. We draw only particles separated from the center of 
the primary by 2 x 10^ cm. The time is indicated at the bottom left in each panel. The time is the initiation time of the carbon detonation 
in the cases of the raw temperature (two left panels) and the smoothed temperature (two right panels). Dashed curves indicate O.IRq 
from the center of the primary. 
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Fig. 17.— Distribution of merger ejecta at t = 500 s in model 5.5M. The ejecta are colored according to their Uej, defined in equation (11). 
All ejecta are projected on the orbital and x-z planes in the top and bottom panels, respectively. 
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Fig. 18.— Distribution of the mass density of particles in model 5.5M at t = 500 s. They are mapped into a one-dimensional profile; 
the horizontal axis indicates the spherical radius. The vertical dotted lines indicate a distance of O.li^Q and l.Oi^Q from the center of the 
merger remnant. 
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Fig. 19.— Mass density and temperature of particles at t = 500 s in model 5.5M. The vertical axes indicate raw and smoothed temperatures 
in the left and right panels, respectively. The dashed curves indicate contours with tec = tu- Above the curves, tec < tu- 
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Fig. 20.— Time evolution of total luminosity of merger shell at the free expansion phase, Sedov phase, and snowplow phase in the case 
of F^k,ej = 3.2 X 10^^ erg. 
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